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ABSTRACT 


Recent  work  on  Ion-electron  recombination  has  shown  that  a  model 

Involving  recombination  Into  highly  excited  bound  levels  through  electron- 

electron-ion  collisions  and  transitions  between  the  bound  levels  through 

electron-atom  collisions  a^  well  as  the  usual  radiative  processes.  Is  In 

good  agreement  with  experiment.  By  making  absolute  intensity  measurements 

of  cll  the  lines  of  helium  emitted  from  the  decaying  plasma  of  an  arc  Jet 

In  the  visible  and  near  ultraviolet,  the  number  densities  of  the  excited 

states  of  helium  have  been  calculated.  When  Interpreted  by  means  of  the 

colllslonal-radlatlve  recombination  model,  these  measurements  give  values 

for  K  .  ,  the  colllslonal  de-excltatlon  rate  constant  of  an  excited 
n,n-l  * 

helium  atom  from  level  n  to  n-1.  The  constancy  of  K  .  over  a 

^  n^n-l 

12  14  3 

range  of  electron  density  from  4  X  10  to  4  X  10  per  cm  indicates 

the  validity  of  the  model,  and  the  values  of  R  ,  are  in  reasonable 

agreement  with  a  recent  theory  of  classical  excitation  cross  sections. 

A  square  root  dependence  of  this  cross  section  on  the  excess  energy  Is 

Indicated,  since  the  values  of  appear  to  be  Independent  of 

temperature.  It  Is  then  found  that  Q2  3  *  ^  ^3  4  ”  X  ^E, 

and  Qy  e  •  380  X  •Je*  where  Q  is  the  averaged  excitation  cross 

,  o  n,  nri 

section  from  the  levels  with  principal  quantum  number  n  to  the  levels 

-15  3 

with  principal  quantum  number  n  -I*  1  In  units  10  cm  ,  and  E  Is 
the  excess  energy  In  ev.  These  cross  sections  have  an  estimated  error 
of  about  307,  In  the  threshold  energy  range  to  about  0.3  ev. 
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INTRODUCTION 

Recent  work  on  electron  recombination  with  monatomic  ions  has  rather 

successfully  explained  the  increase  in  the  recombination  rate  from  the 

values  given  by  radiative  capture  for  low  density  plasmas  to  the  much 

higher  values  generally  observed  in  laboratory  plasmas^  with  densities 

ranging  from  10^^  to  10^^^  ions/cm^. The  two  important  processes 

not  taken  into  account  in  the  older  work  on  the  problem  are  electron 

capture  and  transitions  between  excited  atomic  levels,  both  with  a  free 

electron  providing  the  energy  balance.  A  full  analysis  requires  knowledge 

of  the  cross  sections  for  these  collision  phenomena  for  practically  all  of 

the  excited  levels  of  the  atom,  and  lack  of  this  information  has  been  the 

principal  stumbling  block  in  developing  this  theory.  However,  it  has  been 

found  that  classical  expressions  for  the  ionization  and  excitation  cross 

sections  of  an  atom  by  an  electron  give  recombination  rates  in  quite 

reasonable  agreement  with  experiment.  The  theory  also  predicts  the 

number  of  atoms  in  each  excited  state,  and  one  of  the  principal  verifi* 

(2) 

cations  of  the  theory  rests  on  the  spectroscopic  measurement  of  the 

density  of  these  excited  atoms  In  helium,  and  their  agreement  with  the 
theory.  In  the  present  experiments  on  the  decaying  plasma  of  an  arc  jet 
this  has  been  turned  around,  and  the  spectroscopic  measurements  of  the 
density  of  the  excited  states  of  helium  have  been  used  to  calculate  the 
collisional  transition  rates.  These  are  found  to  be  in  reasonable  agree* 
ment  with  transitional  rates  calculated  using  recently  derived  classical 
cross  sections, in  which  the  motion  of  the  orbital  electron  is  taken 
into  account.  Detailed  calculations  of  the  recombination  rate  using 
these  cross  sections  have  been  given  in  reference  4,  and  some  general 
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simplifications  in  the  theory  are  discussed  in  references  2  and  3«  First 
we  will  discuss  the  plasma  source  and  the  spectroscopic  equipment,  and 
then  the  technique  used  to  obtain  the  density  of  an  excited  atomic  level 
from  the  intensity  of  an  atomic  line.  A  general  discussion  of  these 
experimentally  determined  densities  is  followed  by  a  deduction  of  the 
collisional  de*excitation  rate  constants  for  the  excited  levels,  a 
comparison  of  these  rate  constants  with  those  calculated  using  classical 
cross  sections,  and  finally  the  calculation  of  the  cross  sections  con* 
sistent  with  the  experimental  data. 
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EXPERIMENTAL  EQUIPMENT 

A  schematic  outline  o£  the  plasma  jet  wind  tunnel and  the  spectro¬ 
graph  is  shown  in  Fig.  1.  A  l/4  inch  tungsten  cathode  was  placed  about 
two  inches  axially  from  a  copper  nozzle  which  served  as  the  anode.  The 
nozzle  had  a  throat  diameter  of  0.44  inches^  and  all  parts  were  water 
cooled.  With  a  helium  flow  rate  of  0.05  gm  moles/sec  and  an  arc  current 
of  500  amperes  I  the  pressure  in  the  arc  chamber  was  about  90  mm  Hg  and 
the  arc  voltage  was  about  33  volts.  In  most  of  the  experiments  a  free 
jet  expansion  of  the  helium  plasma  was  used,  and  the  nozzle  ended  at  the 
sonic  throat.  The  main  body  of  the  wind  tunnel,  which  contained  a 
traverse  meclianism  for  positioning  various  probes  in  the  flow  and  a 
quartz  window  for  the  spectrograph,  was  maintained  at  about  0.25  mm  Hg 
by  a  steam  ejector  pumping  system. 

The  problem  of  a  free  jet  expansion  in  an  inviscid  fluid  has  been 
analyzed, and  the  solution  indicates  a  rapid  acceleration  of  the  flow 
to  a  high  Mach  number,  followed  by  a  normal  shock  wave.  The  lower  the 
pressure  in  the  chamber  in  which  the  expansion  occurs,  the  higher  the 
Mach  number  reached  before  the  shock  wave.  Unfortunately,  because  of 
its  small  atomic  mass  the  helium  flow  has  a  low  Reynolds  number,  so  that 
the  viscous  corrections  are  large,  and  in  fact  under  the  conditions  in 
these  experiments  inviscid  flow  theory  cannot  be  applied  to  the  expanding 
flow  to  determine  the  Mach  number  and  other  characteristics.  This  was 
borne  out  by  a  comparison  of  the  axial  variation  of  the  spectroscopically 
measured  electron  densities  with  the  gas  densities  deduced  from  measured 
Impact  pressures  and  isentropic  flow  theory.  This  measurement  indicated 
that  the  Mach  number  never  reached  values  much  higher  than  3. 
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It  was  not  possible  to  run  at  higher  flow  rates  to  reduce  these  vis¬ 
cosity  effects,  as  the  arc  voltage  would  break  lattons  of  as  much 

as  20  volts  amplitude  at  a  frequency  of  a  few  kilocycles^,  accompanied  by 
similar  oscillations  in  the  luminosity  of  the  plasma «  At  low  flow  rates 
all  oscillations  (up  to  about  100  megacycles)  were  either  absent  or  small. 

The  temperature  of  the  helium  at  the  throat  of  the  orifice  can  be 
estimated  from  the  ratio  of  arc  chamber  pressure  with  and  without  arc 
ope rat ion «  It  can  be  shown  that  the  mass  flow  through  a  sonic  throat  is 
proportional  to  p^/  ‘^1  where  p^  and  are  the  stagnation  pressure 

and  temperature 0  and  a  comparison  of  the  measured  stagnation  pressures 
with  and  without  arc  operation  gave  a  stagnation  temperature  of  about 
7000®K.  The  temperature^  density^  and  pressure  in  the  free  jet  expansion 
drop  rapidly 5  although  not  as  rapidly  as  predicted  by  isentropic  flow 
theory 3  until  the  shock  wave  is  reached ^  about  foui  inches  downstream. 

After  the  shock  wava^  which  Is  quite  diffuse^,  the  plasma  had  a  diameter 
of  about  six  inches  and  appeared  quite  uniform  to  the  end  of  the  tunnel ^ 
four  feet  farther  downstream.  When  the  arc  was  operating  '*well*%  this 
plasma  was  remarkably  stationary  and  symmetric^  giving  the  appearance 
of  a  large  glow  discharge  tube.  At  other  times  the  jet  was  less  symmetric^ 
as  the  arc  showed  a  tendency  to  attach  itself  to  one  side  of  the  anode^^ 
and  considerable  flickering  could  be  noticed  in  the  boundaries  of  the  Jet. 

The  spectrograph  used  was  a  Jarrell  Ash  3,4  meter  Ebert  type^  which 
can  take  photographic  spectra^  or  as  was  the  case  in  these  experiments, 
can  be  fitted  with  an  exit  slit  and  a  photomultiplier  and  be  operated 
as  a  direct  recording  instrument.  Wavelength  coverage  was  obtained  by 
rotation  of  the  grating,  with  readout  on  a  counter  directly  in  Angstroms. 
The  radiation  In  a  cross  section  of  the  plasma  about  3.3  mm  high  and 
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0,5  mm  wide  (In  the  direction  of  flow)  was  thus  incident  on  the  slit  of 
the  spectrograph.  The  plane  mirror  swiveled  on  ball  bearings  about  two 
axes  with  calibrated  micrometer  drives^  so  that  any  region  of  the  plasma 
within  about  eight  inches  of  the  orifice  could  be  observed.  This  plane 
mirror  also  rotated  sufficiently  so  that  various  discharge  tubes  and 
radiation  standards  could  be  focused  on  the  spectrograph. 

A  survey  of  the  radiation  from  the  helium  plasma  on  photographic 
plates  revealed  only  neutral  helium  atoms,  with  the  Rydberg  series  merging 
into  a  continuum  at  about  n  «  12  for  the  region  of  the  jet  close  to  the 
orifice.  Farther  downstream  the  Rydberg  series  disappeared  at  about 

4 

n  »  18 and  there  was  no  observable  continuum.  With  long  exposure  (10 
times  that  necessary  to  observe  the  stronger  lines  of  helium),  N2  and 
molecular  bands  were  observed,  principally  in  the  boundary  of  the  jet, 
and  hydrogen  and  argon  i Ines  were  observed  more  uniformly  across  the 
plasma.  The  nitrogen  bands  were  probably  due  to  air  leaks  in  the  tunnel, 
while  the  hydrogen  and  argon  lines  indicated  that  the  gas  purity  was 
better  than  100  parts  per  million.  In  particular,  no  tungsten  or  copper 
from  the  arc  was  found. 

The  argon  impurity  lines  in  helium  show  a  pecularity.  When  the  arc 
is  run  on  pure  argon,  only  the  lines  belonging  to  argon  I,  the  neutral 
atoraj)  are  found.  This  is  presumably  due,  as  in  helium,  to  the  fact  that 
the  more  highly  ionized  species  recombine  very  rapidly,  before  the  plasma 
reaches  the  point  of  obsei'vation.  However,  the  argon  impurity  lines  in 
helium  are  all  of  argon  II  (singly  ionized),  with  no  argon  I  observable. 

It  seems  as  if  some  process  must  very  rapidly  excite  the  ground  state 
argon  ions  before  they  have  a  chance  to  recombine.  One  possibility  is 
that  collisions  of  the  second  kind  with  helium  atoms  in  the  n  «  2  levels 
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may  excite  the  ground  state  argon  ions.  Because  of  resonance  radiation 
trapping,  all  of  the  n  «  2  levels  of  helium  will  have  fairly  high  densities, 
although  considerably  less  than  the  electron  density  (because  of  collisional 
de-excitation) ,  Since  none  of  the  observed  argon  ion  lines  had  an  excita¬ 
tion  energy  greater  than  21.2  volts,  the  highest  energy  of  the  n  2  helium 
levels,  this  process  would  seem  to  be  indicated.  In  order  to  suppress  the 
argon  I  radiation,  it  is  necessary  that  the  rate  of  excitation  of  the  argon 
ions  be  considerably  greater  than  their  rate  of  recombination,  Upon 
estimating  the  recombination  rate  (from  the  theory  presented  in  this  paper) 
and  assuming  the  resonance  radiation  trapping  to  be  complete  so  that  the 
density  of  helium  atoms  in  the  n  »  2  levels  can  be  calculated  assuming 
collisional  de-excitation,  it  is  found  that  the  cross  section  for  excitation 

of  an  argon  ion  in  a  collision  with  an  excited  helium  atom  must  be  consider- 
-14  2 

ably  greater  than  10  cm  .  In  view  of  such  a  large  cross  section,  this 
explanation  must  be  regarded  as  only  tentative c 

The  measurement  of  the  radiation  intensity  from  a  region  of  plasma 
in  the  Jet  cannot  be  made  directly,  since  observations  on  the  centerline 
of  the  plasma  jet,  for  instance,  include  contributions  from  all  regions 
out  to  the  boundary.  If  the  jet  is  assumed  to  be  cylindrically  symmetric, 
the  radiation  reaching  the  spectrograph  is  given  by  the  Abel  integral 
equation,  whose  inversion  is  well  known.  This  inversion  requires  the 
complete  profile  of  the  intensity  of  each  line  across  the  jet,  thus  the 
intensity  of  each  helium  line  was  recorded  as  the  movable  plane  mirror, 
driven  by  a  motor,  performed  a  transverse  scan  of  the  jet.  A  Moseley 
x«y  recorder,  with  x-axis  coupled  to  the  mirror,  was  used  to  record  this 
scan.  The  exit  slit  of  the  spectrometer  was  set  at  0.71  mm  and  the  entrance 
slit  at  0.13  mm,  so  that  quite  uniform  response  was  achieved  over  a 
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wavelength  range  of  3.4  A»  for  reasons  discussed  later.  This  spectral 
response  was  determined  by  scanning  a  singlet  helium  line  from  a  discharge 
tube.  (It  might  be  noted  that  adjustment  of  the  RCA  1F28  photomultiplier 
tube  behind  the  exit  slit  was  required  to  find  a  region  of  the  photo¬ 
cathode  with  uniform  response.) 

The  intensity  calibration  of  the  spectrograph-photomultiplier  system 
was  made  with  a  tungsten  ribbon  filament  lamp,  viewed  with  the  same  opti¬ 
cal  arrangement  as  the  plasma  except  for  the  quartz  window  in  the  tunnel. 
The  primary  standard,  made  by  General  Electric  (#30a/t  24-3),  and  equipped 
with  a  quartz  window,  was  supplied  by  the  Bureau  of  Standards  and  cali¬ 
brated  for  spectral  intensity  from  2500  A  to  7500  A,  when  run  at  35.0 
amperes  AC.  This  lamp  was  used  only  occasionally,  and  another  identical 
lamp  calibrated  against  the  primary  standard  was  used  for  normal  operation. 
The  measurement  of  current  had  to  be  quite  precise,  as  an  error  of  0.3% 
in  current  gave  a  noticeable  change  in  intensity.  The  ammeter  was 
periodically  checked  against  a  laboratory  standard  accurate  to  0.1%. 
Scattered  light  from  the  internal  parts  of  the  spectrograph  necessitated 
the  use  of  a  U.V.  transmitting,  visible  absorbing  filter  for  calibration 
below  3400  A.  With  this  filter,  Corning  7-54,  the  calibration  was 
extended  to  2500  A,  with  about  65%  scattered  light  at  2500  A.  A  red  trans¬ 
mitting  filter.  Corning  2-73,  was  used  to  eliminate  the  second  order 
radiation  for  wavelengths  larger  than  6000  A.  The  whole  system  had  less 
than  2%  variation  in  sensitivity  over  periods  of  several  days.  The  photo¬ 
multiplier  was  powered  by  batteries  and  left  on  continuously,  and  only 
currents  less  than  10  ^  amperes  were  drawn  from  it  to  minimize  fatigue 


effects. 
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A  complete  set  of  measurements  of  all  obsexrvable  lines  of  helium 
required  about  1-1/2  hours ^  and  during  this  time  the  plasma  Jet  was  gen¬ 
erally  quite  steady.  The  intensit>  of  any  given  line  of  helium  varied 
by  less  than  5%  if  the  arc  current  was  maintained  constant,  while  the 
arc  voltage  would  typically  drop  about  10%.  The  largest  errors  in  the 
final  results  probably  arise  from  the  fact  that  small  errors  in  measure¬ 
ment,  or  in  lack  of  symmetry  of  the  jet,  have  a  large  effect  on  the  Abel 
in teg  ra 1  inve  r s ion « 

If  there  is  no  self  absorption  the  equation  describing  the  radiation 
J(x)  observed  by  the  spectrograph  in  a  predetermined  wavelength  interval 
is  given  by 

X  ^  ' 

where  I(r)  is  the  intensity  of  emission  of  the  Jet,  in  watts/cm  ,  assumed 
to  be  a  function  of  radius  only,  and  J(x)  is  the  energy  detected  by  the 
spectrograph  in  watts/cm  ,  when  aimed  at  distance  x  from  the  Jet  center. 
(See  Fig.  Al)  The  inversion  of  this  Abel  equation  is  given  by 

R  dJiiidx 

“  «  J  ^2  _  i  xjl 

r  ^  ' 

(8  9  10) 

Numerical  solutions  based  on  Eq.  (2)'  ^  ^  have  been  described  in  the 

literature.  The  experimental  points  which  represent  the  curve  of  J(x) 
obtained  in  the  transverse  scan  of  the  jet  are  fitted  with  an  inter¬ 
polating  function,  either  directly  or  by  some  least-square  curve  fitting 
technique,  so  as  to  preserve  maximum  accuracy.  Three  different  inter¬ 
polating  methods  for  the  experimental  points  were  tried,  a  linear  fit  to 
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/Q\  /Q\ 

two  points  at  a  time,  ^  a  cubic  fit  to  four  points  at  a  time,  ^  and  a 
six  term  Fourier  cosine  series  fitted  to  twenty  points, By  use  of  a 
function  which  approximated  our  experimental  curves  in  shape,  and  whose 
analytic  inversion  was  obtained,  the  accuracies  of  these  various  methods 
were  tested,  while  varying  the  total  number  of  points  representing  J(x) 
from  5  to  20,  and  the  accuracy  of  each  point  from  four  significant  figures 
to  two  significant  figures.  In  Appendix  A  a  detailed  discussion  of  these 
procedures  is  give:  For  two  significant  figures,  whicn  represent  the 
relative  accuracy  of  our  experimental  measurements,  it  was  found  that  the 
cubic  fit  to  about  five  points  was  as  accurate  as  any  method  tried,  and 
therefore  was  used  for  the  data  reduction  in  these  experiments.  A  relative 
error  of  17  in  the  first  two  or  three  points  near  the  centerline  of  J(x) 
would  result  in  about  a  bX  error  in  I(r)  on  the  centerline.  Since  the 
experimental  curves  were,  at  best,  of  this  order  of  symmetry,  a  error 
is  about  the  best  obtainable  on  the  centerline.  The  two  halves  of  the 
experimental  curves  J(x)  were  averaged  together  before  the  data  reduction 
in  most  of  the  experiments.  The  numerical  work  was  done  on  an  Ijy  !•  "  ) 
computer,  with  Eq.  (2)  put  in  the  form 

N 

I  -  )  T  z  (3) 

m  Nh  ^  n  n,m 

n*0 

where  r  «  mh,  x  »  nh,  and  Z  was,  for  the  cubic  and  linear  fits, 

n,m 

a  triangular  matrix.  For  the  six  term  Fourier  series  expansion  Z 

n,  m 

is  a  20  term  square  matrix. 

The  intensity  of  an  atomic  line,  in  ergs/(cm  -sec-steradian),  is 


given  by 
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where  Is  the  number/cm  of  atoms  In  the  upper  state,  is  the 

Einstein  emission  coefficient  from  state  u  to  state  £  in  sec  ^ 

(averaged  over  all  sub-levels),  and  hv  is  the  energy  of  the  quantum 
of  radiation  in  ergs.  If  the  signal  from  the  photomultiplier  of  the 
spectrometer,  when  set  at  wavelength  X  to  measure  radiation  from  states 
u  -  i  ,  is  called  >  then  the  values  of  to  be  used  in  Eq,  (3) 

will  be  given  by 


J 


ul 

n 


S(\) 


(5) 


where  R^^(X)  is  measured  at  points  whose  distance  from  the  center  of  the 
plasma  is  nh.  As  shown  in  Appendix  B,  W  is  the  spectral  band  pass  (the 
spectrometer  has  approximately  uniform  response  over  a  range  of  wavelength 
W,  and  zero  response  outside  of  this  region),  and  S(X),  the  spectral 
sensitivity  is  given  by 

S(X)  « 

Here  I'*(X)  is  the  spectral  steradiancy  of  the  standard  lamp  and 

is  the  signal  measured  when  the  spectrometer  is  focussed  on  the  standard 

lamp. 

To  find  the  desired  quantity,  ,  in  Eq.  (4),  theoretical  values 
of  the  Einstein  coefficient,  A^^  ,  must  be  used.  Generally,  theoretical 
workers  give  the  oscillator  strengths,  f^^  ,  which  are  related  to  the 
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Einstein  coefficients  by  the  expression 


O  2  2 

8n  e 


me  X 


2  ^iu 


(7) 


where  and  are  the  multiplicities  of  the  upper  and  lower  states, 

and  e,  m,  c,  and  X  are  in  cgs  units.  The  values  found  by  the  coulomb 

f  12) 

expansion  method  of  Bates  and  Damgaard  and  the  variatioual  technique 

(13) 

of  Trefftz,  et.al.'  agreed  quite  closely,  and  were  used  as  the  basis  of 
our  values.  The  tables  of  Bates  and  Damgaard  give  values  for  all  six 
series  of  lines  in  the  visible  up  to  n  «  5  or  6,  and  Trefftz,  et.al. 
calculated  a  number  of  the  astrophysically  important  oscillator  strengths, 
with  some  values  at  n  -  10,  These  values  were  extrapolated  graphically 
to  give  all  oscillator  strengths  up  to  n  «  15  by  using  the  hydrogen 
oscillator  strengths  for  the  same  series  as  a  guide,  and  the  resulting 
values  used  in  this  work  are  given  in  Table  I. 

The  final  expression  for  given  by  combining  Eqs.  (3),  (4), 

(5),  and  (7),  and  results  in  „ 


6.35  X  10  X 


N  /g 
u' 


""  I 

n.^0 


nm 


S(X)  X  EfX)  X  X  h 


(8) 


where  W  is  the  spectral  band  pass  in  Angstroms^  is  the  photo- 

-8 

multiplier  current  in  units  of  10  amp  for  the  transition  u  i,  and 

S(X)  is  the  spectral  sensitivity  of  the  spectrometer  in  units  of 

"8  2 
10  amp/ (microwatts/ [steradian-millimicron-mm  ]).  E  is  given  by 

2 

0.662  X  gj/^  >  with  X  in  microns,  f^^  is  the  oscillator  strength, 
and  h  is  the  interval  in  cm  between  the  values  of  J  • 
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RESULTS 

Measurements  were  made  on  all  the  lines  of  helium  that  could  be 

observed,  and  some  of  the  transitions  responsible  for  these  lines  are 

indicated  on  the  energy  level  diagram  shown  in  Fig.  2.  The  centerline 

values  of  from  Eq.  (8)  were  plotted  on  a  logarithmic  scale 

versus  the  energy  level,  as  shown  in  Figs.  3  and  4.  These  graphs  are 

(2) 

similar  to  those  obtained  by  Hinnov  and  Hirschberg,  and,  as  was  the 

case  with  their  results,  all  the  states  of  higher  excitation  seem  to 

fall  on  a  straight  line,  indicating  a  Boltzmann  distribution  whose 

temperature  is  given  by  the  slope.  The  Boltzmann  distribution  arises 

from  the  fact  that  transitions  between  the  levels  induced  by  electron 

collisions  are  much  more  rapid  than  the  net  rate  of  recombination,  a 

result  leading  to  an  estimate  of  the  cross  section  in  agreement  with 

the  results  to  be  presented  later. 

For  the  lower  levels,  even  though  they  deviate  from  a  Boltzmann 

distribution,  it  appears  that  all  levels  with  the  same  principal  quantum 

number  are  approximately  maintained  with  relative  populations  given  by 

3 

a  temperature  equal  to  that  found  from  the  higher  levels.  The  3  S  and 
3 

4  S  levels  seem  to  deviate  the  most  from  this  relationship,  but  an  error 
in  the  oscillator  strengths  for  these  levels  could  change  the  results. 

It  seems,  then,  that  collisionally  induced  transitions  between  the  triplet 
and  singlet  states  belonging  to  the  same  quantum  numbers  must  be  fairly 
frequent,  compared  to  the  rate  at  which  transitions  occur  from  other 
levels.  The  general  overall  agreement  of  results  from  the  different 
series  for  the  higher  levels  indicates  that  the  oscillator  strengths 
used  in  the  data  reduction  are  accurate  to  at  least  about  10%,  and  that 
the  measurements  themselves  are  accurate  to  about  107., 
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The*  high  c;uantutn  number  levels  are  also  assumed  to  have  rapid  tiai.yi- 
tions  induced  by  electron  collisions  with  levels  lying  above  the  ionization 
limit,  that  is,  with  the  free  electrons.  Thus,  the  whole  system  of  upper 
levels  is  in  equilibrium  with  the  free  electrons,  and  the  Saha  equation. 


N  N, 
e  i 


IlIs. 

8u 
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2nm  kT 
e 


f  ^  -E^kT 
)  ^ 


(9) 


can  be  used  to  give  the  density  of  the  free  electrons.  Upon  setting 
“  ^i  inserting  numbers,  Eq.  (9)  becomes 


(1.2  X  10^^  T 


3/2  ,  1/2 


where  N 

e 


is  in  electrons/cm" 


T  is  in  ev,  and 


-E  /kT 


The  quantity  N J%^  is  found  graphically  by  the  extrapolation  of  the 
upper  levels  in  Figs.  3  and  4  with  a  straight  line  to  the  E  «  0  axis. 

The  electron  densities  and  temperatures  obtained  are  given  on  the  figures. 

The  Abel  integral  inversion,  Eq.  (2),  can  be  considered  as  finding 
the  mean  depth  of  plasma,  which  is  divided  into  the  intensity  measured 
by  the  spectrometer  to  give  the  intensity  per  unit  volume  of  the  plasma. 
For  the  higher  principal  quantum  numbers  this  mean  plasma  depth  should  be 
a  smooth  function  of  quantum  number  only,  since  all  levels  belonging  to  a 
given  quantum  number  are  in  equilibrium  with  each  other.  A  plot  was  made 
of  this  mean  depth  versus  principal  quantum  number,  and  a  smoothing  curve 
passed  through  the  points.  The  results  of  the  Abel  integral  inversion 
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were  then  adjusted  so  that  the  mean  depth  for  each  spectral  line  was  that 
given  by  the  smoothing  curve.  In  some  instances  this  procedure  yielded 
appreciably  less  scatter  in  the  plots  of  vs.  energy.  Figure  3 

was  treated  in  this  manner,  while  Fig.  4  was  not.  In  most  cases  the  points 
were  not  adjusted  more  than  107,,  and  also  the  mean  depth  of  the  plasma  did 
not  generally  vary  by  more  than  307,  over  all  the  lines. 

The  lines  with  high  quantum  number  lying  near  the  continuum  approach 
3  A  in  width,  and  the  assumption  that  the  line  width  is  less  than  the 
spectral  slit  width,  made  in  deriving  Eq,  (5),  is  no  longer  true.  For 
these  lines  a  spectral  scan  was  made,  with  the  spectrometer  focused  on 
the  axis  of  the  plasma  column,  and  the  area  under  the  resulting  curves 
was  measured  with  a  planimeter  to  find  the  intensity.  The  mean  plasma 
depth  as  described  previously  was  then  used  to  find  the  number  of  atoms 
in  the  upper  levels.  These  values  are  the  points  in  Figs.  3  and  4  for 
n  greater  than  10,  and  while  they  seem  to  be  a  little  below  the  equilib¬ 
rium  value,  this  is  probably  due  to  experimental  error. 

Two  independent  verifications  of  the  electron  density  were  attempted. 
First,  the  Stark  broadening  of  the  line  (n  «  4),  when  about  17o 

hydrogen  was  added  to  the  helium,  was  measured.  The  ha If -width  of  the 
3'P  helium  line  was  also  measured  and  compared  with  a  discharge  tube  to 
estimate  the  temperature  of  the  atoms  in  the  plasma,  and  the  half-width 
of  minus  the  instrumental  and  Doppler  half-widths  was  called  the 

half-width  due  to  Stark  broadening.  The  electron  density  was  then  cal¬ 
culated  using  the  results  on  hydrogen  line  broadening  of  Griem,  Kolb, 

Q4) 

and  Shen,  and  Table  II  gives  some  results.  The  rather  poor  agreement 
is  not  unreasonable  in  view  of  the  large  errors  in  the  half-width  measure¬ 
ments,  lack  of  correction  for  spatial  nonuniformity,  and  the  simple 
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subtraction  of  the  various  half -widths .  Particularly  In  the  lower  density 
case,  the  Instrumental  half-width  makes  a  large  correction.  Secondly,  the 
Inglls-Teller  relation^^^^ 


N 
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2a 


0.027 

3  7.5 
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0 


•8 

where  a^  Is  the  Bohr  radius,  0.53  X  10  cm,  and  n  Is  the  quantum 
number  of  the  highest  observable  discrete  level  before  the  lines  merge 
Into  a  continuum,  was  used.  This  formula  was  derived  on  the  basis  of 
the  Stark  broadening  of  the  highly  excited  states  of  hydrogen,  but  should 
be  expected  to  apply  to  helium.  Table  II  shows  the  results  of  this  cal¬ 
culation,  and  It  can  be  seen  that  the  densities  are  In  good  i^greement 
with  the  values  from  the  Saha  equation.  Consideration  of  random  errors 
leads  to  a  probable  error  of  about  157.  In  the  values  of  electron  density 
from  the  Saha  equation. 

The  analysis  of  the  recombination  process  follows  the  same  model  as 
employed  In  references  2,  3,  and  4.  A  rate  equation  Is  written  down  for 
each  level  of  principal  quantum  number  n,  taking  Into  account  both 
radiative  processes  and  colllslonal  processes  Involving  electrons. 

These  equations  are  given  by 
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where  K  is  the  collisional  rate  constant  for  the  transition  from 
m,n 

level  m  to  n.  The  absorption  of  radiation  is  not  explicitly  included^ 
but  will  later  be  taken  into  account  by  dropping  all  radiative  transitions 
to  the  ground  level,  and  the  sums  over  the  discrete  levels  can  be  general¬ 
ized  to  include  the  continuum.  For  the  plasmas  considered  here,  is 

much  larger  than  any  except  the  ground  level,  and  the  parameters 

of  the  plasma  such  as  temperature  and  density  change  slowly  enough  in  time 
so  that  all  the  dN^/dt  except  dN^/dt  in  the  complete  set  of  Eqs,  (10) 
can  be  set  equal  to  zero.  By  then  making  some  reasonable  approximations 

regarding  the  nature  of  K  and  using  the  experimental  measurements  as 

m,  n 

summarized  in  Figs.  3  and  4,  values  of  K  .  can  be  obtained  for  those 

m,  m*  1 

levels  which  are  important  in  determining  the  recombination  rate.  This 

method  depends  on  knowledge  of  the  Einstein  coefficients,  and  knowledge 

of  the  approximate  behavior  of  K 

m,  n 

Recently  Gryzinski^^^  has  calculated  cross  sections  for  excitation 
of  atomic  levels  by  electron  impact  from  classical  considerations  that 
appear  to  be  in  fair  accord  with  experiment.  The  excitation  rate  can  be 
written  in  terms  of  the  cross  section  as 
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n,m 


r  dN  (E) 

J  -3i-  ''<»  « 


where  is  the  excitation  energy  of  the  level,  dN^(E)/dE  is  the 

electron  energy  distribution,  v(E)  is  the  electron  velocity,  and 

Q  (E)  is  the  excitation  cross  sectlo)i  xhe  principle  of  detailed 
n,m  r  r 

balance,  which  states  that  at  equilibrium  the  forward  and  reverse  rates 
of  all  possible  reactions  are  equals  implies  that 
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N  N  K  «  N  N  K 
e  n  n^m  e  m  m,n 


when  all  the  constituents  of  the  reaction  are  at  thermodynamic  equilibrium. 

Using  the  equilibrium  values  of  N  and  N  at  temperature  T  and  the 

n  m 

Boltzmann  distribution  for  dN^(£)/dE  and  v(E),  the  de-excitation  rate 
coefficient^  in  terms  of  the  excitation  cross  section,  can  be  shovm  to  be 


G  =  g  K 
m,n  m,n 
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In  terms  of  this  expression  Eq.  (10),  with  dN^/dt  *  0,  can  be 


rewritten  as 
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Calculations  of  6  using  Gryzinski's  cross  sections  show  that  the 
m,n  ® 

following  approximations, 


G  A 

0.24 

G  1 

m,m-2 

m,m-l 

G  o 

0.10 

G  ^ 

m,m- j 

m,m-l 

Wl,m-1 

- 

0.03 

G  « 

m,m-l 

18. 


are  obeyed  for  all  n  of  Interest.  Using  these  approximations  and  neg¬ 
lecting  ^1^  other  collisional  transitions,  Eq.  (12)  can 

be  written  as 


^n+ljn  ^^n+l,n  ^  ^nfl,n^ 
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The  factors  D  ,  A  ,  and  L  can  be  determined  from  the  experi- 
n+l,n*  n+l,n*  n  ^ 

mentally  measured  densities  of  the  various  levels  for  n  >  3.  The  method 

of  obtaining  G  then  consists  of  estimating  G-  «  (in  some  cases 

nri ,  n  I  ^ 

this  term  is  negligible),  and  solving  Eq.  (13)  for  G^  y  This  value  of 
G^  ^  is  then  used  to  solve  the  equation  for  n  ■  4,  etc.,  until  the 
terms  ^  become  zero  as  the  density  of  level  n  approaches  equilib¬ 

rium  with  the  free  electrons. 

In  the  calculation  of  the  radiative  factor  in  Eq.  (13)  the  terms 

involving  radiation  to  the  ground  level  are  dropped,  as  the  radiation  is 
trapped.  Assuming  Doppler  broadened  lines  (Stark  and  collision  broadening 
are  negligible),  the  spectral  absorption  coefficient  at  the  line  center 

^  “1  4  .  u  (16) 

in  cm  atm  is  given  by 


P(0) 


(14) 


-2  -1 

where  is  the  integrated  absorption  coefficient  in  cm  atm  and 

related  to  the  absorption  oscillator  strength  by^^^^ 


S-  -  2.38  X  10^  f, 
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In  Eq.  (14)  m  Is  the  atomic  mass,  and  on  Is  the  reciprocal  wavelength 
of  the  transition  in  cm  For  the  helium  resonance  transition  l^S  -  2^P 
at  537  A,  -  0  .072.  Assuming  typical  conditions  in  the  plasma  jet, 

T  *  1000®K  and  1.7  X  10^^  atoms  in  the  ground  state,  the  absorption  at 
line  center  given  by  Eq.  (14)  is  found  to  be  50  cm  Thus  the  mean  free 

path  of  resonance  radiation  at  the  line  peak  is  about  0.2  mm,  so  that  to 
a  good  approximation  in  Eq.  (13)  the  resonance  radiation  can  be  neglected 
compared  to  the  other  terms. 

The  results  of  this  calculation  of  G  ,  are  shown  in  Fig.  5,  as 

n,n-l  ®  * 

well  as  the  values  calculated  using  Gryzinski's  cross  sections  with 
averaged  energies  of  the  helium  levels.  The  value  of  2  j  important 
at  the  higher  electron  densities,  was  adjusted  so  as  to  bring  all  the 
points  into  reasonable  agreement.  Agreement  with  the  classical  theory 
is  quite  good;  however,  it  appears  to  give  values  somewhat  too  large  for 
the  levels  with  n  >  5,  and  somewhat  too  small  for  levels  with  n  <  4. 

The  relative  constancy  of  G  «  over  a  factor  of  100  in  electron 
density  can  be  regarded  as  a  verification  of  the  model  of  recombination. 

The  rate  of  recombination,  according  to  the  model,  is  given  by 


dN 
_ e 

dt 

w3 

with  the  radiative  term  being  dominant  in  all  cases  except  those  with  the 
highest  electron  density.  As  stated  previously,  resonance  radiation  to 
the  ground  level  is  neglected  as  it  is  trapped.  In  Table  III  the  recom¬ 
bination  coefficient  Ct  as  defined  by 
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is  given  for  most  of  the  experiments,  using  the  results  of  Eq«  (15). 

Values  of  <X  interpolated  from  the  table  of  Bates^  Kingston,  and 
(4) 

McWhirter'  for  trapped  resonance  radiation,  and  from  the  graph  given 

(2) 

by  Hinnov  and  Hirschberg,  are  shown  in  Table  III  for  comparison. 

These  values  depend  on  the  experimentally  determined  and  N^. 

The  values  of  BateSs  et.al.,  are  about  a  factor  2  greater  than  the 
experimental  values  from  £q.  (15) a  somewhat  surprising  result. 

Their  calculations  were  based  on  the  same  model  of  recombination  used 
in  the  data  reduction  and  used  Gryzinski's  expression  for  the  collision 
cross  sections,  so  it  would  seem  that  the  two  results  should  agree, 
within  random  errors  of  50%  or  so.  The  use  of  Einstein  coefficients 
for  hydrogen  Instead  of  helium,  and  the  differences  in  the  collision 
cross  section  found  by  the  present  experiment,  should  make  only  a  small 
correction. 

The  results  of  Hinnov  and  Htrschberg  are  a  more  approximate  calcu¬ 
lation,  still  based  on  the  same  model  of  recombination,  using  collisional 
cross  sections  estimated  from  the  classical  ionization  cross  section  and 
the  oscillator  strengths  of  the  allowed  optical  transitions.  They  were 
made  for  helium,  and  agree  with  their  experimental  measurements  on  the 
recombining  helium  plasma  in  a  stellarator  with  an  error  of  less  than 
507o.  Their  values  are  also  generally  somewhat  greater  than  those  of  the 
present  experiment. 

The  temperature  dependence  of  the  de-excitation  rate  constanto 

Eq.  (11),  is  related  to  the  form  of  the  excitation  cross  section, 

Q  (E*tE  )•  If  we  assume  that 
n,m  ex 
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then  it  is  easy  to  show  that 


G 

m^n 


An  examination  of  the  experimental  data  plotted  in  Fig.  5  reveals  no 

systematic  variation  of  G  with  electron  temperature,  indicating 

m,  n 

that  k  ■  1/2  for  the  range  of  electron  energy  around  kT.  With  the 
approximation  that  kT  «  ,  the  expressions  given  in  Table  IV  can 

be  calculated  for  the  average  excitation  cross  sections.  This  excitation 
cross  section  must  be  considered  to  be  a  sum  over  all  the  levels  belonging 
to  the  upper  principal  quantum  number,  and  an  average,  according  to  the 
distribution,  over  all  the  levels  belonging  to  the  lower  quantum  number. 
Theoretical  considerations^^^^  show  that  these  cross  sections  have  a 
square  root  dependence  on  E  near  E  >  0,  and  the  results  given  here 
are  consistent  with  this  behavior  up  to  an  energy  of  about  0.3  ev. 
Assuming,  then,  this  dependence  on  E  from  threshold  up  to  an  energy 
of  about  0.3  ev,  the  values  of  the  proportionality  constants  in  Table  IV 
have  an  estimated  error  of  about  307«. 

The  classical  cross  section  derived  by  Gryzinski^^^  applies  to  the 
energy  transfer  between  two  free  electrons.  The  expression  given  by  his 
Eq.  (23),  when  written  in  the  form  a(E^,^,E2)  d(^),  is  the  cross 
section  for  excitation  of  an  electron  with  energy  E^  to  an  energy 
between  E^^  +  ^  and  E^  +  +  d(^)  ,  in  a  collision  with  another 

electron  of  energy  E2«  The  direction  of  the  velocity  of  electron  1  has 
been  assumed  to  be  isotropic  in  space,  and  collisions  with  long  inter¬ 
action  times  have  been  neglected.  In  applying  the  above  results  to 
excitation  of  an  atom,  the  bound  electron  is  assumed  to  be  in  a  circular 
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Bohr  orbit,  from  which  the  kinetic  energy  can  be  found.  The  cross 

section  for  excitation  of  a  level  with  principal  quantum  number  n  to 
level  m  is  thus  given  by 

On.m  -  -('l -“.'2’ 

where  ^  is  the  excitation  energy^  E  -E  ,  of  the  level,  and  E  is 
the  energy  of  the  level  in  the  same  sexics  with  principal  quantum  number 
m  +  1.  This  form  gives  the  threshold  variation  of  the  cross  section  as 

1/2 

(E^-Z^)  '  in  agreement  with  the  quantum  mechanical  considerations. 

Use  of  this  cross  section  to  calculate  results  in  fair  agreement 

with  the  experiment^  as  already  shown  in  Fig.  5,  It  appears^  though,  that 
the  calculated  value  of  Q2  ^  is  too  small  by  a  factor  of  3^  and  for  the 
higher  quantum  numbers  the  calcuiated  value  may  be  somewhat  too  large. 

It  might  be  expected  froir  correspondence  principle  arguments,  that  the 
classical  cross  section  would  be  asymptotic  to  the  actual  values  at  large 
quantum  number.  In  this  case^  however  the  whole  problem  should  be  re¬ 
considered.  It  is  possible  to  put  the  model  of  recombination  in  a 
classical  framework,  and  then  changing  these  equations  to  a  quantum 
mechanical  form  by  use  of  the  correspondence  principle  may  not  lead  to 
exactly  the  same  cross  section  as  indicated  above. 

The  cross  section  for  excitation  of  hydrogen  from  the  ground  IS 

level  to  the  2P  level  has  recently  been  measured  by  atomic  beam 
( 181 

methods'  '  in  the  energy  range  near  thresholdy  and  the  results  fit 
the  curve 


(16) 
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for  energies  up  to  3  ev.  licsults  for  excitation  from  IS  to  2S  have 
also  been  determined,  and  although  they  have  not  been  shown  to  fit 

a  curve  of  the  form  of  Eq.  (16),  the  cross  section  appears  to  be  somewhat 
less  than  that  of  Eq.  (16).  It  is  interesting  to  note  that  this  value  for 
hydrogen  agrees  quite  well  with  Gryzinski’s  theory.  If  it  is  true  that 
the  excitation  cross  sections  for  hydrogen  and  helium  are  about  the  same 
for  the  higher  quantum  numbers,  then  it  appears  that  Gryzinski^s  theory 
gives  about  the  correct  results  for  ^  factor  3  too  low  for  Q2.3» 

and  again  about  the  correct  result  for  Q3.4  ^4-5’ 


24. 


SUMMARY 


Detailed  measurements  of  the  number  density  of  the  excited  levels  of 

helium  In  a  recombining  plasma  have  been  shown  to  be  In  good  agreement  with 

the  col 11s Iona l*radlatlve  model  of  recombination.  These  measurements  led 

to  a  determination  of  the  average  excitation  cross  section  of  several 

levels  of  helium  In  the  threshold  range  to  0.3  ev  excess  energy  which  are 

proportional  to  the  square  root  of  the  excess  energy,  as  suggested  by 

quantum  mechanical  calculations.  These  cross  sections  are  also  In 

approximate  agreement  with  values  calculated  by  the  classical  expression 

of  Gryzinski. The  tables  of  recombination  coefficient  calculated  by 

(4) 

Bates,  Kingston,  and  McWhirter  are  in  approximate  agreement  with 
the  recombination  coefficient  inferred  from  experiment,  as  expected, 
since  Bates,  et.al..  used  Gryzinskl’s  cross  section  in  their  calculations. 
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TABLE  I 

HELIUM  OSCIUATOR  STRENGTHS 
All  values  of  in  units  of  10 


Princ ipal 
Quantum 
Number 

3  3 

2'’p-n'^D 

2^P-n^D 

3  3 

2-^P-n'^S 

2^P-n^S 

23s-n3p 

2^S-n^P 

3 

62.3 

72,5 

5.7} 

16.5 

7.64 

4.64 

A 

12.3 

12.1 

2.08 

5.08 

1.16 

0.824 

5 

4.67 

4.31 

1.19 

2.29 

0.385 

0.300 

6 

2.33 

2.05 

0,709 

1,26 

0.190 

0.151 

7 

1.32 

1,15 

0.414 

0.720 

0,114 

0.087 

8 

0.849 

0.740 

0.268 

0.465 

0.071 

0.055 

9 

0.580 

0.490 

0,186 

0,  311 

0.048 

0.037 

10 

0,408 

0.351 

0.131 

0.223 

0.034 

0.026 

11 

0.300 

0.280 

0.099 

0,168 

0.025 

12 

0.225 

0.195 

0.077 

0.129 

0.018 

13 

0.175 

0.152 

0.060 

0.103 

0.014 

14 

0,140 

0.120 

0.049 

0.083 

15 

0.112 

0.098 

0,039 

0.069 

TABLE  II 
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TABLE  III 

COMPARISON  OF  RECOMBINATION  RATE 


CX  in  10  cm  sec 


N 

cm 

T 

e 

ev 

from  present 
experiment 

from  Bates, 
et.  al. 

from  Hinnov 
and  Hirschberg 

12 

4.2  X  10 

0.15 

5.8 

14 

13 

13 

1.6  X  10^-^ 

0.28 

1.9 

3.3 

1.5 

13 

1.9  X  10 

0.15 

30 

68 

58 

13 

2.6  X  10 

0.13 

89 

150 

130 

13 

6.6  X  10 

0.17 

34 

65 

111 

14 

4.4  X  10 

0.29 

25 

33 

13 

30, 


TABLE  IV 

EXCITATION  CROSS  SECTION  FOR  HELIUM  ATOMS 

Averaged  over  all  levels  with  principal  quantum  number  n 
and  summed  over  all  levels  with  principal  quantum  number  n+1. 

-15  2 

in  units  of  10  cm  with  E  in  ev 

‘=2.3  *  ’-3  X 

S.4  -  '3  X  </e 

Q,  c  (E+E  )  =  380  X  Ve 

^4,5  '  ex' 
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FIG.  4.  DENSITY  OF  EXCITED  STATES  OF  HELIUH 
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F16.  5.  COLLlSIOmiL  DB-EXCITA.TIOM  BATE  CONSTANTS  FOR  HELIUM 


HTD  8001 
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APPENDIX  A 

THE  NUMERICAL  SOLUTION  OF  THE  ABEL  INTEGRAL  EQUATION 


As  shown  in  Fig.  Al,  the  spectrograph  observes  all  the  radiation 
along  line  L  emitted  in  a  solid  angle  Q  .  This  radiation  will  be 
given  by 


J(x) 


r(x,y)  dy 


(Al) 


-y 


O 


where  I(x,y)  is  the  intensity  per  unit  solid  angle  of  emission  at  a 

point  x,y  in  the  plasma,  the  limits  extending  to  the  edge  of  the  plasma. 

3  2 

If  I(x,y)  is  in  watts/cm  stcradian,  then  J(x)  will  be  in  watts/cm 

steradian.  If  the  plasma  is  assumed  to  by  cylindrically  symmetric  about 


X  »  y  *  0,  then  I  is  a  function  only  of 
changing  the  variable  of  integration  in  Eq. 


r  =  (x 
(Al)  to 


2 


+  y 


2, 


1/2 


r,  there 


Upon 

results 


J(x) 


rdr 


2)  172 


(A2) 


where  R  gives  the  radius  of  the  plasma  at  the  edge-  This  is  called  a 

( 21 ) 

Volterra  integral  equation  of  the  first  kind.  Since  the  kernel  has 

a  singularity  at  one  of  the  end  points,  the  standard  techniques  of 

solving  such  equations  cannot  be  used  until  a  cransforma c ion  is  applied 

^22) 

to  Eq.  (A2)  which  removes  this  singularity.  When  this  is  done,  and 

the  inversion  transformed  back  to  the  original  functions,  there  results 


dx 


dx 

7772 


I(r) 


(A3) 
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Unless  J(x)  is  an  appropriate  analytic  functioni  it  is  necessary  to  use 
numerical  techniques  to  solve  this  equation. 

When  J(x)  is  given  as  a  table ,  of  experimentally  determined 

points,  the  most  obvious  technique  is  to  break  Eq.  (A3)  into  a  sum  of 
integrals,  with  each  integral  covering  the  region  between  two  entries  in 
the  table  of  J^.  Some  analytic  interpolating  function  is  then  used  to 
represent  J(x)  in  the  region  between  the  points,  the  integral  in  Eq. 

(A3)  is  performed,  and  the  result  can  be  represented  as 

N 

I  -  e  )  J  Z  ,  (A4) 

m  h  n  nm  ’  ' 

n«0 

where  the  radius  of  the  plasma  has  been  divided  into  N  intervals  of  size 
h,  such  that  x  «  nh  and  r  «  mh. 

Choosing  the  interp>olating  function  as  a  polynomial  of  third  ord^r, 
we  let 

J(x)  ■  Ax^+Bx^  +  Cx  +  D  (A5) 

n  n  n  n 


represent  the  value  of  J(x)  between  x  «  nh  and  x  «  (n^-l)h.  The 

constants  are  evaluated  from  J  J,  J.-,,  except  that  for 

n-1*  n*  n+1'  n+2' 

n  ■  0  the  condition  *  0  at  x  •  0  is  substituted  for  J  , . 

dx  n-1 

Substituting  Eq.  (A5)  in  Eq,  (A3)  there  results 


I(mh) 


1 

7T 


N  (n+l)h 

E  I 


n«0 


[3A  x^  +  2B  X  +  C  ]dx 
n  n  n 


(A6) 


When  the  integrals  in  Eq.  (A6)  are  carried  out,  and  the  values  of 

B  ,  and  C  in  terms  of  J  are  substituted,  it  is  found  that  i...  *  terms 
n’  n  n 

in  J^,  Jj^,  etc.,  can  be  collected.  Eq.  (A6)  can  then  be  put  in  the  form 
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N 

T(tnh)  «  i  )  J  Z  .  (A6) 

n»0 

The  values  of  Z  were  calculated  on  an  IBM  1620.  and  the  values  up 
n ,  m 

to  It  »  20  foi  m  »  0  aie  given  in  Table  AI.  Tliese  values  are  the  same 

(9) 

as  those  given  by  Bockasten,  except  that  he  chose  a  different  manner  of 
fitting  the  interpolating  function  at  the  edge  of  the  plasma,  giving  some¬ 
what  different  values  for  the  last  two  entries  of  Z  .  The  resulting 
difference  in  the  calculated  value  of  I(r)  near  the  centerline  is 
negligible . 

If  a  linear  interpolating  function  is  used  in  place  of  Eq,  (A5),  then 

somewhat  different  values  of  Z  are  found.  Values  for  m  =  0,  taken 

n,m 

from  Nestor  and  Olsen,'"  ^  are  given  in  Table  AI.  If  the  values  of 

are  accurate  enough,  use  of  the  coefficients  Z  derived  from  a  linear 

n ,  m 

rather  than  a  cubic  interpolating  function  will  always  give  less  accurate 

values  of  I  .  However,  as  can  be  seen  from  Table  AI,  A  linear 
- -  m  nra 

introduces  less  cancelation  near  the  center  of  the  values  of  J  than 

n 

Z  cubic,  so  that  if  the  values  of  J  are  not  too  accurate,  and  a 
nm  n 

large  number  of  points  are  used,  more  accurate  results  may  be  obtained  by 
use  of  Z,  .  .  In  this  case,  though,  either  gi*aphical  or  analytic 

smoothing  of  the  points  will  increase  the  accuracy.  A  least  squares 
technique,  applied  to  the  whole  profile  J(x)  as  discussed  by  Freeman 
and  Katz,^^^^  or  to  several  groups  of  points  as  discussed  by  Bockasten, 
can  be  used  to  give  analytic  smoothing.  If  the  values  of  are  suffic¬ 

iently  accurate,  however,  use  of  a  smoothing  technique  may  reduce  the 
accuracy  of  I  . 
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A  smoothing  technique  based  on  the  use  of  an  incomplete  Fourier  cosine 
series  representation  of  was  tried  and  compared  in  accuracy  to 

the  use  of  cubic.  This  method  is  equivalent  to  least  squares  smoothing 

applied  to  the  whole  profile  J(x), 

Equation  (A5)  is  replaced  by 

k 

J(x)  “  ^  ®  R  (A7) 

a»0 

where  f^  are  given  by  the  normal  Fourier  inversion, 

R 

-  I  J(x)  cos  ^  dx  ,  (a  ^  0)  .  (A8) 

o 

Upon  inserting  Eq.  (A7)  in  Eq.  (A3),  there  results, 

k 

°  R  ^  ®  ®a 
a»0 

with 

B^(P)  = 

p  V-  -  ^ 

Equation  (A8)  is  then  evaluated  using  the  trapezoidal  rule  with  the 
given  values  of  J(x)  ■  , 

N 

'a  ■  I  tV^  *  I  Ja  cos  (aycn/N)]. 

n«l 

This  is  substituted  in  Eq.  (A9)  and  the  order  of  summation  reversed, 
resulting  in 


/sin  ajtu 
..2.  .2.1/2 


40. 


with 


and 


I 

m 


1 

h 


I 


J  Z 
n  n,in 


(AlO) 


k 

Z  »  2  )  a  B  (m/N)  cos  (artn/N)  ,  n  /  0 

a»l 


Z 


o,m 


k 

'C - 1 

2^  a  (m/N) 
a*l 


Values  of  Z^  ^  given  by  Shumaker for  k  *  6  and  N  *  20  are 

presented  in  Table  AI,  It  is  seen  that  the  values  of  I  calculated 

m 

with  these  coefficients  are  less  sensitive  to  errors  in  J  . 

n 

The  accuracy  of  the  values  of  calculated  with  these  different 

tables  of  Z^  was  studied  by  trying  them  on  a  function  which  approximated 
the  experimentally  measured  J(x)j  and  whose  inversion  could  be  obtained 
analytically.  The  function  chosen  was 


J(x) 


8  (1  -  x^) 


3/2 


whose  inversion  is 

I(r)  =  6  (1  - 


J(x)  was  divided  into  5^  10,  and  20  points,  and  the  accuracy  was  varied 
by  rounding  off  to  4,  3,  and  2  significant  figures.  For  two  significant 
figures  of  accuraty,  which  was  estimated  to  be  the  error  in  reading  the 
traces  in  these  experiments,  calculation  of  using  20  points  and 

Z^  Fourier  was  comparable  in  accuracy  to  using  five  points  and 
cubic.  The  error  was  about  TU  near  the  centerline  in  this  case.  If  20 
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points  and  cubic  were  used,  the  error  was  much  larger,  about  207e» 

However,  when  given  to  four  significant  figures,  the  error 

in  I  using  20  points  and  Z  Fourier  was  still  about  27.,  while  the 
m  ntn 

error  using  Z  cubic  with  20  points  was  less  than  0.1%.  In  all  cases 
nm 

use  of  Z  linear  resulted  in  less  accuracy  than  use  of  Z  cubic, 
nm  •'ntn 

Some  of  the  experimental  results  were  evaluated  with  Z^^  Fourier 

as  well  as  Z  cubic,  and  no  significant  increase  in  accuracy  was  noted, 
nm 

Probably  the  smoothing  introduced  by  eye  when  reading  the  values  of 
from  the  charts,  coupled  with  division  of  J(x)  into  about  seven  points, 
is  sufficient  to  give  about  maximum  attainable  accuracy. 
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TABLE  AI 


COEFFICIENTS  FOR 

THE  ABEL  INTEGRAL 

INVERSION 

in 

-2 

units  of  10 

n 

Z  ,  cubic 
no* 

interpolation 

Z  ,  linear 
no 

interpolation 

Z  ,  Fourier 
no' 

interpolation 

0 

76.258 

63.661 

8.240 

1 

-58,007 

-42.443 

12.481 

2 

-  5.847 

-  8.488 

3.020 

3 

-  3,395 

-  3.635 

-  6.090 

4 

-  1.970 

-  2.018 

-  9.772 

5 

-  1.269 

-  1.286 

-  6.994 

6 

-  0.883 

-  0,891 

-  1.013 

7 

-  0.649 

-  0.656 

3.403 

8 

"  0.497 

-  0.499 

3.521 

9 

-  0,393 

"  0.395 

0.281 

10 

-  0.318 

-  0.318 

-  2.988 

11 

-  0.263 

-  0.263 

-  3.504 

12 

-  0,222 

-  0.222 

-  1.156 

13 

-  0.188 

-  0.188 

1.748 

14 

-  0,162 

-  0.162 

2.677 

15 

-  0.141 

-  0.141 

]  .026 

16 

-  0.125 

-  0.125 

-  1.548 

17 

-  0.110 

-  0.110 

-  2.749 

18 

-  0.098 

-  0.098 

-  1.695 

19 

-  0.088 

-  0.088 

0.398 

20 

-  0.081 

-  0.081 

... 

43. 


I 


FIG.  Al.  DIAGRAM  OF  PORTION  OF  CYLIHDRICALLY  SYMMETRIC  PLASMA 

VIEWED  BY  SPECTROMETER 


HYD  8002 
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APPENDIX  B 

THE  MEASUREMENT  OF  THE  INTENSITY  OF  AN  ATOMIC  LINE 

Einstein  assumed  that  an  atom  in  an  excited  state  u  has  a  certain 

probability  per  unit  time,  A^^  ,  of  making  a  spontaneous  transition  with 

(23) 

emission  of  radiation  to  state  £  of  lower  energy.  With  the  develop¬ 

ment  of  quantum  theory  it  became  possible  to  compute  from  first  principles 
the  value  of  A^^  for  simple  systems.  The  intensity  of  an  atomic  line  is 
given  by 


u£ 


4n 


- 


(Bl) 


3  3 

in  ergs/(cm  -sec-steradian) ,  with  the  number/cm  of  atoms  in  the 

upper  state  and  E^  -  E^  the  energy  of  the  quantum  of  radiation. 

Absorption  and  induced  emission  are  neglected  in  this  formula.  Thus, 

if  A  «  is  known,  measurement  of  I  .  will  lead  to  a  value  of  N  , 

First  it  is  noted  that  the  wavelength  of  the  transition  from  state 

u  to  state  i  is  not  perfectly  monochromatic.,  but  is  spread  over  a 

wavelength  range  AX  by  various  line  broadening  phenomena.  In  these 

experiments  Doppler  broadening  is  the  most  important  for  the  low  lying 

levels  of  helium,  while  the  electric  microfields  due  to  nearby  electrons 

and  ions  cause  a  large  broadening  of  the  levels  of  helium  near  the 

ionization  limit,  The.'ie  line  broadening  phenomena  do  not  change  the 

mean  radiative  lifetime  of  the  upper  state  appreciably,  so  that  Eq.  (Bl) 

still  holds,  with  I  ,  considered  to  be  all  the  radiation  associated 
*  ui 

with  the  line.  The  spectrometer  responds  to  the  radiation  falling  on 

2 

the  entrance  slit,  J(X)  in  ergs/(cm  -sec-Angstrom-steradian) ,  which 


1 


must  be  first  integrated  over  the  line  width  A\  corresponding  to  tiie 
transition  u  /, 
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K(.  “  / 

and  then  the  Abel  integral  inversion,  Eq.  (2),  must  be  performed  on 
to  find  of  Eq,  (Bl).  The  Abel  integral  inversion  in  effect  finds 

the  mean  depth  D  of  che  emitting  plasma,  which  gives  as 


I 


ujf. 


Ki/'> 


(B3) 


The  magnitude  of  the  photomultiplier  current  when  the  spectrometer 
is  set  at  wavelength  \  can  be  given  by 

R(^)  *  J  8(>--^')  .  (B4) 

where  g(A“\*)  is  the  response  function  of  the  spectrometer  and  is 
different  from  zero  for  only  a  small  range  of  wavelength  around  Its 

shape  is  determined  by  the  entrance  and  exit  slit  widths  and  the  aberra¬ 
tions  of  the  spectrometer,  and  its  magnitude  principally  by  the 
sensitivity  of  the  photomultiplier  and  transmissivity  of  the  spectrometer* 
When  the  spectrometer  is  focused  on  the  standard  tungsten  filament  lamp, 
its  signal  will  be 

R*(>^)  -  J  -  J«(\)  J  g(X-\')  d\'  , 

where  J®(X*)  is  the  known  spectral  steradiancy  of  the  standard  lamp. 
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and  can  be  taken  out  of  the  Integral  since  it  will  be  constant  over  the 
wavelength  range  where  g(\*^c')  differs  from  zero.  Thus  the  sensitivity 
of  the  spec t rome ter ^ 

S(\)  -  R'’(X)/j“(\)  ,  (B6) 


is  determined,  and  a  typical  graph  of  this  versus  wavelength  is  shown  in 
Fig.  Bl.  When  the  spectrometer  is  focused  on  the  plasma,  and  the  wave* 
length  is  set  to  receive  radiation  from  transition  ui,  the  signal  is 


AX 


g(X^-X') 


I 

AX 


d\’ 


(B7) 


where  the  response  function  g(X^-X*)  has  been  assumed  to  have  a  constant 
value  over  the  width  AX*  of  the  line.  This  was  achieved  in  the  present 
experiment  by  having  an  entrance  slit  v;idth  of  0.71  mm,  resulting  in  a 
reasonably  uniform  response  function  over  a  wavelength  range  of  3  A.  If 
in  Eq.  (B7)  the  line  width  is  small  compared  to  changes  in  g(X-X*),  the 
line  acts  as  a  delta  function  in  picking  out  a  value  of  X*,  and  R^|(X) 
will  be  a  faithful  representation  of  g(X*X*).  Upon  producing  a  curve  of 
R^^(X)  by  scanning  the  spectrometer  slowly  over  a  narrow  line,  and 
integrating  the  resulting  trace  with  a  planimeter,  there  results 


f  ^  y  g(X-X')  dX' 


R  -(X  ) 
u/'  o' 


(B8) 


where  value  of  the  curve  at  its  flat  top,  W  Is  called 

the  spectral  slit  width,  in  Angstroms,  and  is  independent  of  X  within 
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the  experimental  error.  It  can  now  be  seen  that 

P  R  .(\  )  X  W 

■  J  ■  -it--  ■  <”> 

AX  ° 

which  is  identical  with  Eq.  (5). 

W  was  determined  using  the  singlet  lines  of  helium  from  a  discharge 
tube  in  several  regions  of  the  spectrum^  and  did  not  vary  by  more  than 
5%,  The  width  of  these  lines  is  about  0,07  A,  while  g(X-\*)  has  a 
flat  top  about  3  A  wide,  so  that  the  approximations  leading  to  Eq.  (B8) 
are  well  satisfied.  Equation  (B9)  will  correctly  describe  the  intensity 
of  a  line  emitted  from  the  plasma  as  long  as  the  width  of  the  line  is 
somewhat  less  than  3  A. 
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